Microscopic theory of weak pseudogap behavior in the underdoped cuprate 
superconductors I: General theory and quasiparticle properties 



OO 
On 
On 

u 



> 

On 
(N 

^- 

o 

00 
ON 

a 

C 

o 
o 



13 



Jorg Schmalian, a David Pines, a,b and Branko Stojkovic fc 
a University of Illinois at Urbana- Champaign, Loomis Laboratory of Physics, 1110 W. Green, Urbana, IL, 61801 
h Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, NM, 87545 

(February 1, 2008) 

We use a novel solution of the spin fermion model which is valid in the quasi-static limit nT 3> w s f , 
found in the intermediate (pseudoscaling) regime of the magnetic phase diagram of cuprate super- 
conductors, to obtain results for the temperature and doping dependence of the single particle 
spectral density, the electron-spin fluctuation vertex function, and the low frequency dynamical spin 
susceptibility. The resulting strong anisotropy of the spectral density and the vertex function lead 
to the qualitatively different behavior of hot (around k = (7T, 0)) and cold (around k = (n /2,n/2)) 
quasiparticles seen in ARPES experiments. We find that the broad high energy features found in 
ARPES measurements of the spectral density of the underdoped cuprate superconductors are de- 
termined by strong antiferromagnetic (AF) correlations and incoherent precursor effects of an SDW 
state, with reduced renormalized effective coupling constant. Due to this transfer of spectral weight 
to higher energies, the low frequency spectral weight of hot states is strongly reduced but couples 
very strongly to the spin excitations of the system. For realistic values of the antiferromagnetic cor- 
relation length, their Fermi surface changes its general shape only slightly but the strong scattering 
of hot states makes the Fermi surface crossing invisible above a pseudogap temperature T*. The 
electron spin-fluctuation vertex function, i.e. the effective interaction of low energy quasiparticles 
and spin degrees of freedom, is found to be strongly anisotropic and enhanced for hot quasiparticles; 
the corresponding charge-fluctuation vertex is considerably diminished. We thus demonstrate that, 
once established, strong AF correlations act to reduce substantially the effective electron-phonon 
coupling constant in cuprate superconductors. 
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I. INTRODUCTION 



In addition to their high transition temperatures and 
the d x 2_ y 2 symmetry of their superconducting state, 
the cuprate superconductors possess a remarkable range 
of normal state anomalies. Seen first as charge re- 
sponse anomalies in transport, Raman, and optical ex- 
periments, and subsequently as spin response anoma- 
lies in nuclear magnetic resonance (NMR) and inelas- 
tic neutron scattering (INS) experiments, recent specific 
heat and angular resolved photoemission spectroscopy 
(ARPES) experiments have shown that these anoma- 
lies are accompanied by, and may indeed originate in, 
anomalous planar quasiparticle behavior. It is conve- 
nient to discuss the temperature and doping dependence 
of this "uniformly" anomalous behavior in terms of the 
schematic phase diagram shown in Fig. 1. There one sees 
that overdoped and underdoped systems may be distin- 
guished by the extent to which these exhibit crossover 
behavior in the normal state: underdoped systems ex- 
hibit two distinct crossovers in normal state behavior 
before going superconducting, while overdoped systems 
pass directly from a single class of anomalous normal 
state behavior to the superconducting state M. In a 
broader perspective of this schematic phase diagram, 
which is applicable to the YBa 2 Cu307_,5, YBa 2 Cu40s, 
Bi 2 Sr 2 CaiCu20 8 +5, HgBa 2 Cu0 4+ 5, HgBa 2 Ca2Cu 3 8 



and Tl 2 Ba2Ca2Cu30io systems, if one defines an opti- 
mally doped system as that which possesses the highest 
superconducting transition temperature within a given 
family, then optimally doped systems are in fact under- 
doped. 

Attempts to understand the different regimes of this 
phase diagram have been based on strong magnetic pre- 
cursors [0-R] , the formation of dynamical charge modula- 



tions in form of stripes |11| 
Cooper pairs above T c 



, the appearance of preformed 
10[ , or the separation of spin 



and charge degrees of freedom [|12|.[13| . 

A phase diagram similar to Fig. 1 was independently 
derived from studies of the charge response by Hwang, 
Batlogg, and their collaborators jlij and from an analysis 
of the low frequency NMR experiments |l5f| by Barzykin 
and Pines Jig] , The latter authors identified the upper 
crossover temperature, T cr , from measurements of the 
uniform susceptibility, \ , in Knight shift experiments, 
which show that for underdoped systems \ possesses a 
maximum at a temperature T cr , which in underdoped 
systems, increases rapidly from T c as the doping level is 
reduced. The fall-off in susceptibility for temperatures 
below T cr was first studied in detail by Alloul et al. [[l7| , 
and led Friedel |p|] to propose that it might arise from a 
near spin density wave (SDW) instability; he coined the 
term pseudogap to explain its behavior, in analogy to 
the quasiparticle pseudogap seen in charge density wave 
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(CDW) systems. Barzykin and Pines identified further 
crossover behavior in this pseudogap regime by examin- 
ing the behavior of the 63 Cu nuclear spin-lattice relax- 
ation time, 63 Xi , as the temperature was reduced below 
T cr . They noted that between T cr and a lower crossover 
temperature, T*, the product, 63 XiT decreases linearly 
in temperature, while shortly below T* this product has 
a minimum, followed by an increase as the temperature 
is further lowered, an increase which is strongly sug- 
gestive of gap-like behavior. They proposed that these 
two crossovers were accompanied by changes in dynam- 
ical scaling behavior which could be measured directly 
if NMR measurements of 63 Ti could be accompanied by 
measurements of the spin-echo decay time, 63 T2G. Above 
T cr they argued that the ratio, 63 TiT/ 63 T 2 2 G , would be 
independent of temperature, a result equivalent to ar- 
guing that the characteristic energy of the spin fluctua- 
tions, oj s { , would be proportional to the inverse square 
of the antiferromagnetic correlation length, £. Between 
T cr and T* they proposed that the ratio, 63 TiT/ 63 T 2G 
would be independent of temperature, which means that 
an underdoped system would exhibit z = 1 scaling be- 
havior, i.e. u> s { would be proportional to ; below T* 
they found that the increase in w s f would be accompanied 
by a freezing out of the temperature-dependent antifer- 
romagnetic correlations; i.e. £ , which was proportional 
to a + bT between T cr and T*, would approach a con- 
stant. This behavior has recently been confirmed in NMR 
measurements on YBa2Cu40s by Curro et al. Jl9) ] while 
2 = 1 pseudoscaling behavior has been found in INS ex- 
periments on La1.86Sro.14Cu O4 by Aeppli et al. [ pp[ . 

Because pseudogap behavior is found both between 
T cr and T*, and between T* and T c , the terms weak 
pseudogap and strong pseudogap behavior were coined 
to distinguish between the two regimes plf . Thus in the 
weak pseudogap regime one finds z = 1 pseudoscaling 
(because the scaling behavior is not universal) behavior, 
with both uj s f and f -1 exhibiting linear in T behavior, 
while the rapid increase in 63 TiT, or what is equivalent, 
u>si, found below T* suggests that strong pseudogap is an 
appropriate descriptor for this behavior. 

An alternative perspective on weak and strong pseudo- 
gap behavior comes from ARPES |^,^3| and tunneling 
experiments [ p4[ , which focus directly on single parti- 
cle excitations. Above T*, ARPES experiments show 
that the spectral density of quasiparticles located near 
the (71", 0) part of the Brillouin zone, develops a high en- 
ergy feature, a result which suggests that the transfer 
of spectral weight from low energies to high energies for 
part of the quasiparticle spectrum may be the physical 
origin of the weak pseudogap behavior seen in NMR ex- 
periments. Below T*, ARPES experiments disclose the 
presence of a leading-edge gap, a momentum-dependent 
shift of the lowest binding energy relative to the chemical 
potential by an amount up to 30 meV for quasiparticles 
near (ir, 0); it seems natural to associate the strong pseu- 



dogap behavior seen in the NMR experiments with this 
leading edge gap. Recent tunneling experiments have 
shown that both the high energy feature (i.e. weak pseu- 
dogap behavior) and the strong single particle pseudogap 
can also be observed in the tunneling conductance, with 
the high energy feature occurring primarily in the occu- 
pied part of the spectrum. 

Strong pseudogap behavior is also seen in specific heat, 
d.c. transport, optical experiments, and Raman experi- 
ments. Below T*, a reduced scattering rate for frequen- 
cies w < 7rT* has been extracted from the optical con- 
ductivity using a single band picture (25| . This suggests 
that excitations in the pseudogap regime are more coher- 
ent then expected by extrapolation from higher temper- 
atures. This point of view is supported by recent Raman 
experiments p^ , p7| which observe in the Bi g -channel, 
sensitive to single particle states around k = (n, 0), a 
suppression of the broad incoherent Raman continuum 
and a rather sharp structure at about twice the single 
particle gap of ARPES experiments Q . Resistivity mea- 
surements also show that below T* the systems gets more 
conducting than one would have expected from the linear 
resistivity at higher temperatures 

It is natural to believe that the pseudogap in the spin 
damping, as observed in 63 TiT measurements, the sin- 
gle particle pseudogap of ARPES and tunneling exper- 
iments, and the pseudogap of the scattering rate are 
closely related and must be understood simultaneously. 
Furthermore, it is essential for any theory of the strong 
pseudogap to account properly for the already existent 
anomalies above T*, because they are likely caused by 
the same underlying effective interactions. As can be 
seen by inspection of Fig. [l], strong pseudogap behavior 
in underdoped cuprates only occurs once the system has 
passed the weak pseudogap state. 

In this paper we concentrate on the weak pseudogap 
(pseudoscaling) regime above and give a more detailed 
account of the preliminary results we obtained using a 
spin-fluctuation model of normal state behavior [El]. It 
is our aim to provide a quantitative understanding for 
T > of quantities reflecting strong pseudogap behav- 
ior below We derive a novel solution of the spin 
fermion model in the quasi-static limit irT 3> w s f, rele- 
vant for the intermediate weak pseudogap regime of the 
phase diagram. We demonstrate that the broad high 
energy features of the spectral density found in ARPES 
measurements of underdoped cuprates are determined by 
strong antiferromagnetic correlations and incoherent pre- 
cursor effects of an SDW state. The spectral density at 
the Fermi energy and the electron-spin fluctuation ver- 
tex function are strongly anisotropic, leading to qualita- 
tively different behavior of hot (around k = (tt, 0)) and 
cold (around k = (tt/2, 7t/2)) momentum states, whereas 
the Fermi surface itself changes only slightly. We present 
new results for the effective interaction of quasiparticles 
with spin and charge collective modes. In distinction to 
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the strong coupling of hot quasiparticles to spin excita- 
tions, we demonstrate that their renormalized coupling to 
charge degrees of freedom, including phonons, is strongly 
suppressed. Finally, we show that the onset temperature, 
T cr , of weak pseudogap behavior is determined by the 
strength, £, of the AF correlations. 

Our theory also allows us to investigate the low fre- 
quency spin and charge response functions. In a sub- 
sequent publication, we will discuss the suppression of 
the spin damping and further generic changes in low fre- 
quency magnetic behavior seen in NMR experiments as 
well as the optical response, particularly as far as the 
Bi ff -Raman continuum is concerned. 

The paper is organized as follows. In the next sec- 
tion we summarize important findings of ARPES results 
which will later on be explained by our theory of the 
weak pseudogap regime. Next, we give the basic concept 
of the spin fluctuation model and derive the spin fermion 
model. In the following, fourth, section we discuss in de- 
tail our solution of the spin fermion model in the quasi- 
static limit, with particular attention to the new physics 
of the spin fermion model for intermediate coupling. Our 
solution is obtained by the complete summation of the 
perturbation series, and is motivated in part by a theory 
for one dimensional charge density wave systems devel- 
oped by Sadovskii p9[ . We have extended his theory 
to the case of two spatial dimensions and isotropic spin 
fluctuations and, in so doing, found that we could avoid 
several technical problems of the earlier approach. Tech- 
nical details of the rules we used for computing diagrams 
are presented in Appendix |A| and B. Readers not inter- 
ested in these technical aspects can skip the theory sec- 
tion and should be able to follow the discussion of our 
results for the spectral density and vertex functions in 
the fifth and sixth sections, respectively. In particular, 
results for the single particle properties are discussed at 
length and compared with ARPES experiments. Finally 
our theory for the weak pseudogap regime is summarized 
in the last section, where we also consider the physics of 
the strong pseudogap state and summarize some predic- 
tions and consequences of our theory. We argue that a 
proper description of the higher temperature weak pseu- 
dogap regime is essential for a further investigation of the 
low temperature strong pseudogap state and argue that 
the strong pseudogap state and precursors in the pairing 
channel are the quantum manifestation of strong anti- 
ferromagnetic correlations whereas the spin density wave 
precursors are the classical manifestation of it. 



II. ARPES EXPERIMENTS 

ARPES experiments offer a powerful probe of the 
quasiparticle properties of cuprates. Since they provide 
unusually strong experimental constraints for any theory 



of optimally doped and underdoped cuprate supercon- 
ductors, we summarize in this section the main experi- 
mental results obtained by this experimental technique 

In Fig. ||, we show ARPES spectra close to the mo- 
mentum k = (-7T, 0), for two different doping concentra- 
tions While for the overdoped, T c = 78 K, sample a 
rather sharp peak occurs, which crosses the Fermi energy, 
the spectral density of the underdoped, T c = 88 K, sam- 
ple exhibits instead a very broad maximum at approx- 
imately 200 meV. Thus, the entire line-shape changes 
character as the doping is reduced. The other important 
difference between the two charge carrier concentrations 
is the appearance of the leading edge gap (LEG), i.e. a 
shift of the lowest binding energy relative to the chemical 
potential, for the underdoped system. This LEG varies 
between 20 and 30 meV and is therefore hardly visible in 
Fig. H , but is discussed in detail in Ref. |p2| , p3| . 

In addition to this strong doping dependence, the 
spectral function of underdoped systems is also very 
anisotropic in momentum space, as can be seen in Fig. [|. 
Here, the position of local maxima of the spectral func- 
tion along certain high symmetry lines of the Brillouin 
zone is shown for an overdoped and underdoped sys- 
tem. This is usually done because the maxima of the 
spectral density correspond to the position of the quasi- 
particle energy. However, as we discuss in detail below, 
this interpretation is not correct in underdoped systems 
for momentum states close to (tt, 0) where the line-shape 
changes qualitatively Close to (ir, 0), one sees for over- 
doped systems, in agreement with Fig. 0, a peak at low 
binding energy, which crosses the Fermi energy between 
(ir, 0) and (ir, 7r), whereas the 200 meV high energy fea- 
ture is the only visible structure for the underdoped sys- 
tem. It is flat and seems even repelled from the Fermi 
energy between (tt,0) and (jr, it). The situation is dif- 
ferent for momentum states along the diagonal, where a 
rather sharp peak crosses the Fermi energy between (0, 0) 
and (7r/2,7r/2); the velocity of the latter states, seen in 
the slope of their dispersion in Fig. pi is independent of 
the doping value, while no LEG has been observed for 
those quasiparticle states. 

In Ref. H|, the authors constructed the Fermi sur- 
face for the two doping regimes by determining the k- 
points where a maximum of the spectral function crosses 
the Fermi energy. Their results are re-plotted in Fig. ^. 
Consistent with Fig. || a large Fermi surface occurs for 
the overdoped material, whereas only a small Fermi sur- 
face sector close to the diagonal could be identified in 
the underdoped case. Even though this appears to be 
in agreement with the formation of a hole pocket closed 
around (ir/2,ir /2), with reduced intensity on the other 
half of the pocket, ARPES data below the superconduct- 
ing transition temperature, shown in Fig. ||[ show that for 
momenta close to (tt, 0), a sharp peak appears at lower 
binding energy. This behavior, for the underdoped case 
is completely consistent with a large Fermi surface which 
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is only gapped due to the superconducting state. The ob- 
vious question arises: how could a transformation from 
a small to a large Fermi surface occur on entering the 
superconducting state? 

We should also mention that in Ref. |30|, the authors 
showed that the two energy scales (the LEG and the high 
energy feature) behave as function of doping in a fashion 
which is quite reminiscent of the two temperature scales, 

and T cr , shown in Fig. El This leads us immediately to 
two conjectures: 1. There is a relationship of the physics 
of the upper crossover temperature T cr and the high en- 
ergy feature, as well as between the strong pseudogap 
temperature T» and the LEG. 2. As a strong pseudogap 
state is impossible without a weak pseudogap state at 
higher temperatures, the LEG can only appear after the 
system has established the high energy features. 

As noted above, these fascinating experimental results 
represent a set of very strong constraints for the micro- 
scopic description of underdoped cuprates we develop be- 
low. 



III. THE SPIN FLUCTUATION MODEL 

The nearly antiferromagnetic Fermi liquid (NAFL) 
model ||^,|3j| of the cuprates offers a possible explana- 
tion for the observed weak and strong pseudogap behav- 
ior. It is based on the spin fluctuation model, in which 
the magnetic interaction between the quasiparticlcs of 
the Cu02 planes is responsible for the anomalous normal 
state properties and the superconducting state with high 
T c and d x 2_ y 2 pairing state ]3^ , |3^ |. In a recent letter, 
we have shown how the weak pseudogap regime can be 
understood within this NAFL-scenario [0 . 

In common with many other approaches, within the 
spin fluctuation model the planar quasiparticles are as- 
sumed to be characterized by a starting spectrum which 
reflects their barely itinerant character, and which takes 
into account both nearest neighbor and next nearest 
neighbor hopping, according to 

£k = — 2i(cos k x + cos k y ) — At' cos k x cos k y — fj, , (1) 

where t, the nearest neighbor hopping term, ~ 0.25 
eV, while the next nearest neighbor hopping term, 
may vary between t' w — 0.45t for YBa2Cu306+,5 and 
t' fa -0.25< for La 2 _ x Sr 2; Cu04. 

In distinction to many other models, the spin fluc- 
tuation model starts from the ansatz that the highly 
anisotropic effective planar quasiparticlc interaction mir- 
rors the dynamical spin susceptibility |3q| , 



V. 



NAFL 
cff 



(q,w) = ,9 2 X q (o;) 



(3) 



XqH ~l + £ 2 (q-Q) 2 -^' 
peaked near Q = (ir.Tr), via: 



(2) 



an ansatz which enables us to construct directly a theory 
which focuses solely on the relevant low energy degrees 
of freedom. In Eq. |3|, g is the coupling constant char- 
acterizing the interaction strength of the planar quasi- 
particles with their own collective spin excitations. In 
this model, changes in quasiparticle behavior both re- 
flect and bring about the measured changes in spin dy- 
namics. The dynamic susceptibility, Eq. ||, was intro- 
duced by Millis, Monien, and Pines |35| to explain NMR 
experiments, which can be used to determine the cor- 
relation length, £, the constant scale factor, a, and the 
energy scale u> s t, which characterizes the over-damped 
nature of the spin excitations. It follows from the ex- 
perimental data that the static staggered spin suscepti- 
bility xq = a £ 2 is large compared to the uniform spin 
susceptibility, xo, and the relaxational mode energy cor- 
respondingly small compared to the planar quasiparticle 
band width Jl6|,[56| . For optimally doped and underdoped 
systems one finds that over a considerable regime of tem- 
peratures, 



(4) 



and it is only as T falls below T* that u> s { becomes 
comparable to and eventually larger than ttT. In de- 
tail, between T, and T cr one finds: uj si /(irT) w 0.17 
for YBa 2 Cu 4 8 and 0J ai /(irT) w 0.14 for YBa 2 Cu 3 6 . 6 3 
rather independent of T |l6| . As a result of this compar- 
atively low characteristic energy found in the weak pseu- 
dogap region, the spin system, for q ~ Q, is thermally 
excited and behaves quasi-statically |2l| ; the quasiparti- 
cles see a spin system which acts like a static deformation 
potential, a behavior which is no longer found below T* 
where ui s f increases rapidly |]l6| and the lowest energy 
scale is the temperature itself. 

Since the dynamical spin susceptibility x q (w) peaks at 
wave vectors close to (ir, tt), two different kinds of quasi- 
particles emerge 37,3a]: hot quasiparticles with 



(5) 



located close to those momentum points on the Fermi 
surface which can be connected by Q, feel the full ef- 
fects of the interaction of Eq. (Ep ; cold quasiparticles with 
|£k — ^k+cjl > v/£, located not far from the diagonals, 
\k x \ — \k y \, feel a "normal" interaction. In Fig. |^, we 
show the Fermi surface in the first quarter of the BZ and 
indicate the evolution with £ of its hot regions, which ful- 
fill Eq. ||, by a thick line. Note that even for a correlation 
length £ = 1 a different behavior along the diagonal and 
away from it is expected. For larger values of £, the hot 
regions become smaller while their effective interaction 
increases. Close to T c , typical values for £ of underdoped 
but superconducting cuprates are 2 < £ < 8, depend- 
ing on doping concentration (Tq |; v is the magnitude of 
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a typical Fermi velocity in the corresponding momentum 
regions. 

The distinct lifetimes of hot and cold quasiparticles 
can be obtained from transport experiments: a detailed 
analysis shows that, due to the almost singular inter- 
action, the behavior of the hot quasiparticles is highly 
anomalous, while cold quasiparticles may be character- 
ized as a strongly coupled Landau Fermi Liquid ]3§| ] . The 
presence of incommensurate peaks in the spin fluctuation 
spectrum [^0|,[39| , and hence in the NAFL interaction, al- 
though difficult to calculate, may be expected to amplify 
the role played by hot quasiparticles in the determination 
of system behavior. 

In the spin fluctuation model the anomalous behavior 
of the cuprates is assumed to originate in a strong interac- 
tion between fermionic spins s q = \ X) k(TCT < c k+q<j CT o-o-' c k<T' 
which brings about intermediate range (£ > 1) antiferro- 
magnetic spin correlations and over-damped spin modes. 
Here, the operator c k(7 creates a quasiparticle which con- 
sists of hybridized copper "&(S. x i_yi and oxygen 2p x (y) 
states The quantity of central physical interest is 

the dynamical spin susceptibility 



(6) 



which after Fourier transformation in frequency space 
and analytical continuation to the real axis is assumed 
to take the form, Eq. ^. The intermediate and low en- 
ergy degrees of freedom are characterized by an effective 
action [B2| 



S=- 



P rP 
dr I dr' 



k,<r 



collective-mode interaction, between these two regimes, 
is beyond the scope of this paper. 

The quantities of primary interest to us arc the 
single particle Green's function Gk, ff (r — t') = 
— (Tcu,cr(T)c^ (r')) which provides information about the 
quasiparticle spectral density determined in angular re- 
solved photo-emission experiments, the dynamical spin 
susceptibility itself, and the corresponding charge re- 
sponse functions. As noted above, in calculating these 
quantities for intermediate correlation lengths the inter- 
action between the collective spin modes is irrelevant. In 
appendix [a| we show that under these circumstances the 



Green's function 



Gka(r-T') = (G kikCTCT (r,T'|S) 



(8) 



can be expressed as a Gaussian average (• • •)„ of electron 
propagators with a time dependent magnetic potential 
-^S q (r), with respect to collective bosonic spin 1 vari- 
ables S q (r). The corresponding model is often referred 
to as the spin fermion model. We give in appendix ^ 
the diagrammatic rules of this problem, which will be es- 
sential for the solution of the spin fermion model in the 
quasistatic limit. In the next two sections, we derive new 
expressions for the single particle Green's function, and 
the spin-fermion and charge-fermion vertex functions of 
the quasistatic two dimensional spin fermion model, valid 
for intermediate values of the spin fermion coupling, by 
extending an earlier study by Sadovskii p9fl for one di- 
mensional charge density wave systems. 



IV. THEORY OF THE QUASISTATIC LIMIT 



^3 z2Xq(.T-T')Bq{T) -B-^t') 



(7) 



where G~ k (r — t') = — (<9 T + e^)8{t — r') is the inverse 
of the unperturbed single particle Green's function with 
the bare dispersion, Eq. [j]. In using Eq. ^, we implicitly 
assume that the effect of all other high energy degrees of 
freedom, which are integrated out to obtain the action 
S, do not affect the Fermi liquid character of the quasi- 
particles. In Eq. (t| the effective spin-spin interaction 
is assumed to be fully renormalized; thus it reflects the 
changes in quasiparticle behavior it brings about, and 
can be taken from fits to NMR and INS experiments. 
We will also assume that the spin degrees of freedom 
are completely isotropic and that all three components 
of the spin vector are equally active. In the case of a 
intermediate correlation lengths 1 < £ < 8, this is the 
appropriate description of the spin degrees of freedom. 
Only for much larger £ tx exp(consi./T), does one enter 
the regime in which even without long range order only 
two transverse spin degrees of freedom are active [ ffT| . 
The physics of the crossover, driven by a collective-mode 



We begin this section by first motivating the qua- 
sistatic limit and discussing its physical consequences by 
investigating the second order diagram with respect to 
the coupling constant g. We then present a solution of 
the spin fermion model which is not restricted to the 
weak coupling regime and provides new insight into the 
intermediate coupling behavior relevant for underdoped 
cuprates. 



A. The second order diagram and the static limit 

In second order perturbation theory, the quasiparticle 
self energy is given by: 



l 



lUJ n + W m - £ k +q 



(9) 



If, for a given temperature T, the characteristic frequency 
of the spin excitations w s f is small compared to the in- 
trinsic thermal broadening of the electronic states, the 
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energy transfer iv m of this state due to an inelastic scat- 
tering process is negligible. Furthermore, in the limit 
ttT >• oj s f, Xq(wm) is dominated by the Matsubara fre- 
quency v rn — 0, so Eq. ^ takes the form: 

SkW = ~g 2 Ts( q )-. , (10) 



q 



£k+Q+q 



with g = g aT and 
S{q) 



Physically, this use of a static approximation for the spin 
degrees of freedom reflects the fact that since the fre- 
quency variation of Xq( w ) takes place on the scale w s f, 
once ttT 3> 0J s f, ah relevant collective spin degrees of 
freedom are thermally excited and the phase space re- 
strictions for scattering phenomena due to the quantum 
mechanical nature of the spins are irrelevant. It follows 
that we can then neglect the w-variation of x q (w). 

For the system we study, experiment shows that the 
dominant momentum transfer q of the spin fluctuations 
is close to the antiferromagnetic wave vector Q = (tt, 7r), 
so that we can expand the energy dispersion as 



£k+Q+q ~ £k+Q + V k+ Q • q 



(12) 



with velocity wg+q = cfek+Q/ 'dk a . Note that in distinc- 
tion to a one dimensional problem, the linearization of 
the electron spectrum in two dimensions is not straight- 
forward. In Eq. [l^, we have linearized with respect to the 
transferred momentum q ss Q = (n, tt), an approxima- 
tion which is justified provided q deviates only slightly 
from the antiferromagnetic wave vector Q, i.e. for sys- 
tems with a sufficiently large antiferromagnetic correla- 
tion length £. Therefore, technically is considered to 
be a small quantity and all related momentum integrals 
are evaluated accordingly. On comparing this approx- 
imate treatment with a complete numerical evaluation, 
we find that it can be applied once £ > 1. At k = (n, 0), 
the velocity v£ vanishes and one must take higher order 
terms in q — Q into account. We assume that the physics 
of this van Hove singularity is irrelevant (due to three di- 
mensional effects and the presence of possible additional 
scattering mechanisms) and introduce a lower velocity 
cut off v c ~ (wk)FS- The remaining momentum inte- 
gration can then easily be carried out. It follows, after 
analytical continuation iuj n — > u> + i0 + , that: 



-A 2 



f.TT 

V 2 



(13) 



-arctanh 



w — £k+Q 



VK+q/0 2 + (w - £k+q) 2 



where A 2 = g 2 aT log(l + (£A) 2 ) and A w it is the upper 
cut off of the momentum summation. Since we are tech- 
nically at high temperatures, our results depend on this 



cut off, which is undesirable. We avoid this problem by 
expressing any cut-off dependence of the theory in terms 
of measurable quantities: Thus on using the local mo- 
ment sum rule (S 2 ) = 3T^ m q x q (w m ), we find that A 
can also be expressed as 



A 2 



,9 2 (S 2 )/3. 



(14) 



We therefore can use this expression for A and deter- 
mine (S 2 ) from the experimentally determined suscep- 
tibility Xq(^m) of Eq. |. This guarantees a reasonable 
estimate for the total spectral weight of the spin excita- 
tion spectrum for the spin fluctuation induced scattering 
processes. 

Consider a given k-point on the Fermi surface (£k = 0) . 
If the Fermi surface is such that the momentum trans- 
fer by Q takes you to another Fermi surface point, i.e. 
£k+Q = 0, it follows from Eq. |l3| that for this momen- 
tum state, a so called hot spot, = £k+Q = 0, the real 
part of the self energy decreases like log(u/)/cJ if u> > 
close to the 1/uj behavior which is a signature of precur- 
sor effects of an spin density wave ||. More generally, 
anomalous scattering processes will continue to modify 
the single particle spectrum dramatically for those mo- 
mentum states for which 



kk-£k+q| < v/£. 



(15) 



This entire region of the BZ behaves in qualitatively dif- 
ferent fashion from the rest of the system; it corresponds 
to the definition of hot quasiparticlcs discussed recently 
by Stojkovic and Pines Q. 

We call attention to the fact that only for the hot quasi- 
particles can we justify neglecting the higher Matsubara 
frequencies. For cold quasiparticlcs with |ek — £ k+Q| > 
v/£ the characteristic energy scale of the spin fluctua- 
tions is no longer ui a { but turns out to be ~ w s f£ 2 [p8[ , 
a quantity which is not, in general, small compared to 
7rT. As a result, our approach, while properly account- 
ing for the anomalously large scattering rate and related 
new physics of the hot quasiparticles, will tend to over- 
estimate the scattering rate for cold quasiparticlcs. Put 
another way, differences in behavior between hot and cold 
quasiparticles will be underestimated in our theory. 

In order to make explicit the role played by the pres- 
ence of SDW precursors in the quasistatic regime, we 
evaluate the above momentum integrals within the ap- 
proximation 



S(q) 



r 2 + ? t 



(16) 



where is the projection of q parallel (perpendicular) 
to the velocity Vk+Q. We then obtain: 



E k (w) 



A 2 



:k+Q +ivk +Q /£ 



(17) 



G 



an expression which, apart from a logarithm, has the 
same anomalous behavior as Eq. |l3]. In the limit £ — > oo 
A is the spin density wave gap and the poles of the re- 
sulting Green's function are the two branches of the mean 
field SDW state discussed by Kampf and Schrieffer . 

For the investigation of higher order diagrams in the 
next paragraph, it will be helpful to introduce the follow- 
ing representation of the second order self energy: 

E k (u)=-if f die^x+^Vk+qCt) (18) 
Jo 



where 



Evaluation of the momentum summation yields for A 



^k+Q(i) = 271-1^0 (*«k+q/£) 



(20) 



where Kq is the modified Bessel function. Using the ap- 



proximation of Eq. 16, this simplifies to 



$k+cj(*) 



(21) 



The tendency towards SDW behavior in the quasi- 
static regime so far relies on the applicability of the sec- 
ond order perturbation theory: visible effects can only 
occur once the correlation length exceeds the electronic 
length scale £ G = v/A « 2v/g. In a weak coupling treat- 
ment, the above discussion is applicable only for large 
correlation length: one therefore has to go beyond sec- 
ond order perturbation theory to be certain whether or 
not SDW precursors are relevant for cuprates with inter- 
mediate correlation length. This is possible only if £ G is 
only a few lattice constants; it implies that we have to in- 
vestigate an intermediate coupling regime. Therefore, we 
present in the next paragraph a procedure which enables 
us to sum the entire perturbation series. 



B. Diagram summation in the quasistatic limit 

To evaluate all higher order self energy diagrams within 
the quasistatic limit, we first derive a compact expression 
for an arbitrary diagram and then, as a second step, sum 
all diagrams of the perturbation series to obtain the self 
energy and single particle Green's function. This sum- 
mation is made possible by the fact that many diagrams 
with rather different topology are, apart from a factor 
which describes multiplicity and sign, identical. 

As first shown by Elyutin in the context of optical re- 
sponse in a random radiation field , diagrams can be 
characterized by the sequence of integer numbers {rij}, 
where rij is the number of interaction lines above the j- 
th Green's function; for an example, see Fig. [?]. In the 



following we prove that in the quasistatic regime, dia- 
grams with the same sequence {rij} are proportional to 
each other. The proportionality factor will be determined 
below. 

An arbitrary diagram of order 2N can, up to a con- 
stant, be expressed as: 



~2N 



■S{q 



N 



qi-qjv 
2N-1 



(22) 



where q a = q a — Q and the ((2N — 1) x N) matrix Rj a 
determines whether q Q (a = 1 . . . N) occurs as a momen- 
tum transfer in the j-th Greens function (j — 1 . . . 2N—1) 
of the diagram, i.e. Rj, a = 1 or 0. In general, each di- 
agram is fully characterized by Rj, a - It is important to 
notice that rij is given by the expression, 



N 



(23) 



a=l 



Since each of the momenta q Q of Eq. |22j is separately 
constrained to lie in a region close to q Q ~ Q, we can 
expand: 



iV 



+v k +jQ E fy* (<J« ~ Q) ' ( 24 ) 

a=l 

where we have used the fact that rij is even (odd) if j 
is even (odd) since at each vertex rij changes by ±1 and 
ni = 1. Shifting all momenta q Q — Q — > q Q and intro- 
ducing, as we have done for the second order diagram, 
2N — 1 auxiliary time variables tj, it follows that: 



S^>(k,c) = Hr-V w £ S( qi )...S(q N ) 

qi-qjv 

TT / dtj e rij r~ £k+]Q_Vkj L a =i i? 3-° q 

3=1 



(25) 



with Vk,j = Vk+jQ. In this proper time representation of 
the self energy, the different momentum integrals decou- 
ple; on using Eq. O it follows that 



,(2N) 



(k, W ) = (-i) 



2JV-1 r00 
2 V-l ~2N J~J 



J'=l 

Y\_ i>k+jQ(Rj,atj) 



(I 



A' 



dtj e i ( UJ ~ £k+jQ " >tj 



(26) 



Q = l 



In the last step we used the fact that the momentum 
transfer is sufficiently close to Q that we can neglect con- 
tributions of order (q— Q) 4 in 5(q), i.e. we have assumed 
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that ipk(x + y) oc ipw(x)ipk{y) as follows from the use of 
S(q) in Eq. [if]. This is consistent with the restriction to 
momentum transfers close to q ~ Q which motivated the 
linearization of the electron spectrum in Eq. Eq. |26| 
is only valid for hot spots where the velocities and 
V k+Q are almost perpendicular to each other. For cold 
quasiparticles this condition is not fulfilled, so that our 
theory can only give a qualitative account for their scat- 
tering processes. It has been recently pointed out by 
Tchernyshyov [Q that Eq. ^6] is in fact not fulfilled in 
the original one dimensional solution of Ref . j29| . There- 
fore, the ideas developed in Ref. |29| seem to be much 
more appropriate for our two dimensional case. 



bosons, this problem has been solved by Sadovskii p9[ . 
The generalization (see Appendix |A]) to the case of spin 
fluctuations is not straightforward, because of the addi- 
tional (— l) c factor of crossed spin conserving and spin 
flip lines. 

In appendix |^ we derive the multiplicity of a given 
class of diagrams, i.e. the number of identical diagrams 
of a given order of the perturbation series, by solving the 
problem in the special case £ — > oo and using the fact 
that the combinatorics of the diagrams does not depend 
on this limit. Having determined these multiplicities, it 
is possible to sum the entire perturbation series analyt- 
ically. We find the following recursion relation for the 



Eq. 



Since Rj <a is either or 1, it follows immediately from Green's function G^(u) = G^~°\u): 



N 



3) ' 



(27) 



a = l 



with rij given by Eq. [2^. Inserting this result and col- 
lecting all the prefactors, if follows p4|: 



2N-1 roo 

S( 2iV )(k,w) = H) 2JV - 1 A 2JV H / dt 

.7=1 Jo 



(28) 



which proves that a given diagram of order 2N is fully 
determined by the sequence {rij} as well as provides an 
explicit expression for these diagrams. 

On making use of the simplified evaluation of the mo- 
mentum integrals, Eq. |l6|, it follows with the help of 
Eq. M, that 



£< 2JV )(k, w ) 



2N-1 



,2JV 



II — 77 > ( 29 ) 



a result which is useful in determining the multiplicity of 
a given diagram. 

For the actual evaluation of all diagrams, it is essential 
that for each sequence {rij}, there is a unique mapping 
to a diagram without crossing interaction lines, since for 
each {rij} there exists one and only one diagram with- 
out crossing interaction lines. Note, unique is meant in 
the sense of the topology of a diagram, not whether it 
contains longitudinal or transverse spin fluctuations; for 
details see Appendix B. This is illustrated for two cases in 
Fig. ||, where we show two self energy diagrams of order 
<? 4 and g s , which are, within the quasistatic approxima- 
tion identical apart from a proportionality factor. ^From 
these considerations, it follows that it suffices to sum up 
only the non-crossing diagrams taking into account the 
identical crossing diagrams by their proper multiplicity 
factors. The remaining problem is to determine, for a 
given order in the coupling constant, how many identical 
diagrams exist. In the case of charged and uncharged 



<#Vr = 3k V)" 1 - ^ + iA 2 G^" +1 V) 



O'+i), 



(30) 



with Kj = (J + 2)/3 if j is odd and Kj — j/3 if j is even 



9 k j) M = ~i / dt e^-^+iQ) i> k+jQ (ty . (31) 



and 



Eq. |30| is one of the central results of our theory. This 
recursion relation, closed by G^'{uS) — g^\uj) for some 
large value of L, enables us to calculate the single par- 
ticle spectral function A(k, ui) to arbitrary order 2L in 
the coupling constant g (we use L ~ 10 2 — 10 3 ; Eq. |30| 
converges for L ~ 10 2 ). 



C. Spin susceptibility and vertex function 

Within the quasistatic limit of the effective low en- 
ergy quasiparticlc interaction, we can obtain an exact 
expression for the irreducible part of the dynamical spin 
susceptibility and the electron spin fluctuation vertex. 
Note that we are not able to calculate the total suscep- 
tibility. Since we are assuming that the interaction line 
is given by the fully renormalized spin susceptibility, a 
direct approach would lead to an over counting of dia- 
grams. Therefore, we only calculate the irreducible part 
X^(ioj n ) of the total susceptibility x^fi^n)- The latter 
can be expressed as 



1 - /qXq 



(32) 



where the restoring force, / q , is determined by the renor- 
malization of the spin exchange fermion-fermion interac- 
tion through high energy excitations in all other channels. 
/ q is then related in a non-trivial fashion to the under- 
lying microscopic Hamiltonian of the system and has to 
be considered as an additional input quantity. 

By following a procedure analogous to that in deter- 
mining the Green's function in Eq. pi one can show that 
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the irreducible part of the dynamical spin susceptibility 
is given by 

(iu n ) = (ilf (to;„|S)> o> (33) 

where IIq^(ia; n |S) is the irreducible particle hole propa- 
gator for a given spin field configuration. 

n^(*w w |S) = ~ Tr{a a G k+q , k ,^ m ,„,W|S) 

kk'mm' 

xcr /3 G k /_ q ^ k (i^ m /^,zw m |S)| . (34) 

Here Tr . . . refers only to the trace in spin space and 
= u rn + uj n This result is obtained by neglecting all 
reducible contributions in taking the functional deriva- 
tive with respect to an external time dependent magnetic 
field coupled to the electron spins Sj(r). 

The diagrammatic rules described in appendix |A| for 
the single particle Green's function G k (tj) can be ex- 
tended in straightforward fashion to the spin suscepti- 
bility, which can be expressed in terms of G k (aj) and the 
electron spin fluctuation vertex function: 

Xq("^m) = -5 ^ G k (iw„)G k+q (iw„ + iv m ) 
" k,n 

r k,k+q(^«,iw„-l-zi/ m ). (35) 

Thus, a knowledge of the vertex function gives immediate 
information about the irreducible part of the dynamical 
spin susceptibility. A similar relation exists for the cor- 
responding charge susceptibility. In this paragraph we 
outline the exact determination of r^ k+ (cj + i0 + ,o; + 
v + iO + ) obtained after analytical continuation to the real 
axis. For the determination of the susceptibility on the 
real axis we will also need the analytical continuation 
k+q (w — «0 + , oj + v + i0 + ) which has to be determined 
independently but can be obtained in a similar way. 

As was the case for the electronic Green's function, 
the vertex function is obtained in two steps: first, based 
on purely diagrammatic arguments we obtain a general 
expression for the vertex function in terms of the previ- 
ously determined Green's function and some combinato- 
rial prefactors which take the proper multiplicity of the 
diagrams into account; second, these prefactors are deter- 
mined in the limit £ — > 00. This is possible because the 
combinatorics of the diagrams does not depend on the 
actual value of £. Finally, we obtain a closed expression 
valid for all values of £. 

In the case of the Green's function each diagram was 
proportional to a rainbow diagram. The correspond- 
ing conclusion for the vertex function is that each ver- 
tex diagram is identical to a diagram of the ladder 
approximation and the entire perturbation series can 
be obtained by summing the ladder series with appro- 
priate weighting factors. The proof of this statement 



is almost identical to the corresponding proof for the 
Green's function. Because an arbitrary diagram can be 
related to a ladder diagram, it follows that the vertex 

r £k + <>+*o + ^+^o + ) = r { ^(cu+iO+^+v+iO+) 

can be generated by two renormalized Green's functions 
and an effective vertex r^ 1 k +q (uj+iO + , uj+v+iO + ) which 
includes all those processes not taken into account by one 
spin fluctuation propagator crossing the external bosonic 
line. For the spin vertex, we find the recursion relation, 

r k°k+ q (^ + *0 + , lo + v + z0+) = 1 - n A 2 G^ 1} («) 
4+,( w + ^) r £k+c> + *0+,c + v + z0+) . (36) 

Here, the Green's function G^ (u>), Eq. ||(], takes into ac- 
count that for the diagram under consideration one has 
at least one interaction line above each fermionic propa- 
gator . A comparison with perturbation theory shows 
that the multiplicity factor which enters Eq. [36] is given 
by ri = I . The minus sign in Eq. |36| results from the di- 
agrammatic rules of appendix |A|. Since the higher order 
vertex function can be determined in exactly the same 
way as Eq. [36], one obtains the recursion relation 

r k,k+ q ( w + i0+ ' bJ + v + l0+ ) = 1 - n+iA 2 G£ +1) («) 
xG£J(w + "jrgS'V + *0+, u + v + z0+) , (37) 

which can be evaluated using the Green's functions from 
Eq. |(] and a starting value = 1. In Fig. § the dia- 
grammatic motivation for this recursion relation is given: 
there one sees that the problem is similar to the summa- 
tion of the ladder series for the vertex, with the difference 
that all non-ladder diagrams are taken into account by 
the corresponding weighting factors r; +1 Once these pref- 
actors are known, the vertex function can be determined 
up to arbitrary order of the coupling constant. The n are 
defined diagrammatically by the fact that 3 L Yif=i r i ^ s 
the number of skeleton diagrams of order g 2L which con- 
tribute to the vertex function. [Note, that non-skeleton 
diagrams are diagrams with interaction lines which only 
renormalize the Green's functions.] 

The combinatorial determination of the r; is somewhat 
cumbersome. We proceed by using the general expres- 
sion, Eq.[3?], to calculate the irreducible susceptibility 
given in Eq. |35|, while determining the irreducible sus- 
ceptibility independently in the limit £ — > 00 analytically 
by evaluating the path integral of Eq. On comparing 
these two results for £ — > 00 order by order in the cou- 
pling constant we are able to determine the prefactors r\ . 
On carrying out this calculation for arbitrary momentum 
q, we find r; = I if I is even and r/ = (I + 2)/9 if I is odd. 

This completes the specification of the vertex function, 
Eq.|3^, of the spin fermion model and enables us to calcu- 
late both the irreducible spin susceptibility, Xq(u>), and 
the effective spin fluctuation induced pairing interaction. 
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An identical procedure can be performed for the charge 
vertex IT k+q (lj + i0 + , lj + v + i0 + ) . One obtains in place 
of Eq|37| the result: 



k,k+q 



(lj + i0+, lj + v + j0+) = 1 + n l+1 A 2 G^ +1) (lj) 



xG£J(w + *)rgJi> + l0+ ^ + " + * 0+ ) . ( 38 ) 

since for the charge vertex ri is replaced by — Ki, with m 
as given in the line below Eq. BO. 



V. QUASIPARTICLE PROPERTIES: THEORY 
COMPARED WITH EXPERIMENT 

We consider first the frequency and momentum depen- 
dence of the spectral density, A^lj) = — jlmGk(w), for 
a typical underdoped system. In Fig. [to] we show, in the 
inset, the Fermi surface, defined by those k-points which 
fulfill 



e k + Re£ k (u;) 



(39) 



for lj — 0, for interacting quasiparticles whose bare spec- 
trum is specified by t = — 0.25eV, t' — — 0.35i, at a 
hole doping concentration, = 0.16. In this and all 
subsequent plots we assume g = 0.8 eV, in agreement 
with transport measurements |38|. This corresponds to 
an intermediate regime for the coupling constant since 
it is similar to the total bandwidth. The calculation is 
carried out at a temperature such that £ = 3, which, as 
we shall show, lies in the weak pseudogap regime well 
below T CI . In the main part of Fig. |l^, we show our 
results for Ak(cj)f(u>), where f(u>) is the Fermi function 
for several points on the Fermi surface. A\^(lj) f (lj) , the 
quantity measured in ARPES experiments, is strongly 
anisotropic. For a representative cold quasiparticle (a), 
located close to the diagonal, with \sk — £k+Q| > 
the peak in the spectral density crosses the Fermi sur- 
face. For these quasiparticles, the quasistatic magnetic 
correlations simply produce a thermal broadening of the 
spectrum, as is characteristic of a Landau Fermi liquid 
at small but finite T. 

The situation is completely different for the hot quasi- 
particles at (d) which are located close to (tt, Q). Here, 
|£k — £k+Q| <v/£. A large amount of the spectral weight 
is shifted to higher energies, a shift which gives rise to 
weak pseudogap behavior. As will be discussed below, 
the position of the maximum of this broad feature, which 
represents the incoherent part of the single particle spec- 
trum (i.e., does not correspond to a solution of Eq. 39), 
is similar to the quasiparticle bands of a mean field spin 
density wave state. Thus, even though incoherent in na- 
ture and considerably broadened, this high energy feature 
is the precursor effects of a spin density wave state. A 
second interesting aspect of the calculated hot quasipar- 
ticle spectral density is that although there exists a solu- 
tion of Eq. (|39]) at lj — those quasiparticles [and quite 



generally those near (tt, 0)] do not possess a peak. This 
part of the FS is therefore not observable in an ARPES 
experiment. Experimentally, a FS crossing can only be 
determined if a local maximum of the spectral density 
crosses the Fermi energy. The calculated visible part of 
the FS, where our calculated spectral function exhibits 
a maximum at lj = 0, is shown in Fig. |ll| (thick lines). 
It is in agreement with experiment. While this behavior 
appears to be similar to that expected for a hole pocket, 
below we discuss the important differences between our 
results and a hole pocket scenario. 

The reason for the "disappearance" of pieces of the 
Fermi surface in the weak pseudogap regime is the follow- 
ing. The finite imaginary part of the self energy at lj = 
invalidates, as always for T ^ 0, a rigorous quasiparticle 
picture and can even affect the occurrence of a maximum 
of the spectral density in the solution of Eq. [59|. This is 
what happens for hot quasiparticles in the weak pseudo- 
gap regime. Due to their strong magnetic interaction the 
related large scattering rate causes the hot quasiparticle 
peak to be invisible in the weak pseudogap regime and 
care must be taken to properly interpret the calculated 
Fermi surface. 

Consider now the evolution of the Fermi surface with 
temperature, or what is equivalent, with £. As can be 
seen in Fig. [ll], for £ = 1, the FS is basically unaffected 
by the correlations, a situation very similar to the one ob- 
tained within a self consistent one loop calculation. This 
confirms the result obtained by Monthoux, p§] that ver- 
tex corrections, neglected in the one loop framework, are 
indeed of minor importance for small correlation lengths. 
[Note that while our calculations are based on the fact 
that the dominant momentum transfer occurs near the 
antiferromagnetic wave vector, which implies at least an 
intermediate correlation length £, it is useful to consider 
the limiting case, £ ~ 1, even though in this regime dif- 
ferent theoretical approaches may be turn out to be more 
appropriate.] On increasing £ to values which are realis- 
tic for underdoped cuprates (£ = 2 • • • 8) , we find slight 
changes of the FS-shape for momenta close to (tt, 0) and 
(0, tt); however, the general shape (large FS closed around 
(tt, tt) and equivalent points) remains the same. If one 
further increases £ to values larger then 10 lattice con- 
stants, serious modifications of the FS, caused by a short 
range order induced flattening of the dispersion of the 
quasiparticle solution, begin to occur. This follows from 
the solution of Eq. [39] for finite lj. It is only for such large 
correlation lengths that a hole pocket starts to form along 
the diagonal. Eventually, at some large, but finite value 
of £, our solution gives a closed hole pocket. We conclude 
that for underdoped but still superconducting cuprates 
(with £ < 8), the shape of the FS remains basically un- 
changed, while our theory can potentially describe the 
transition from a large Fermi surface to a situation with 
a hole pocket around (tt/2, tt/2), which may be the case 
very close to the half filling. 



10 



The above results provide a natural explanation for 
what is seen at temperatures below the superconduct- 
ing transition in ARPES experiments on the underdoped 
cuprates: the sudden appearance of a peak in the spec- 
trum of quasiparticles located near (ir,0). According to 
our results, this is to be expected, since as T falls below 
T c , the scattering rate of the hot quasiparticles drops dra- 
matically; the superconducting gap has suppressed the 
strong low frequency scattering processes which rendered 
invisible the peak in the normal state, and a quasiparticlc 
peak emerges. Since this sudden appearance of the quasi- 
particle peak below T c is inexplicable in a hole pocket 
scenario, the ARPES experimental results support the 
large Fermi surface scenario we have set forth above. 

Another interesting aspect of the calculated results 
shown in Fig. [To] is the sudden transition between hot 
and cold quasiparticles, justifying the usefulness of this 
terminology a posteriori. To demonstrate explicitly the 
anisotropy of the spectral function for low frequencies, 
we show, in Fig. [l^, Ak(u> = 0) along the Fermi surface 
as function of the angle 0p = axct&n(ky/k x ) between 
k and the k x axis. Even though no gap occurs in the 
hot quasiparticle spectral density in the weak pseudogap 
regime, the low frequency spectral density is consider- 
ably reduced. It is therefore not possible to consider the 
behavior above the strong pseudogap crossover temper- 
ature, T*, where our theory should apply, as being con- 
ventional. 

We compare, in Fig. [l^, the calculated variation of 
the maximum of A^(cj) in momentum space with the 
ARPES results of Marshall et al. |H| for two different 
doping concentrations. For an overdoped system, we as- 
sumed a correlation length £ = 1 and a charge carrier 
concentration = 0.22. The resulting dispersion cor- 
responds to that of the original tight binding band with 
slightly reduced band width. The plotted maxima for 
£ = 1 all correspond to broadened coherent quasiparticle 
states. We chose t! = — 0.35i leading to a Fermi sur- 
face crossing along the diagonal as well as between (ir, 0) 
and (tt, tt) in agreement with experiments The situation 
is different for an underdoped system, which we assumed 
to have a charge carrier concentration rih — 0.16 and a 
correlation length £ = 3, similar to other underdoped 
but superconducting cuprates. We use the same value 
t' = — 0.35£ for the next nearest neighbor hopping inte- 
gral. Along the diagonal, we still find a Fermi surface 
crossing and, in agreement with experiment, no doping 
dependence of the Fermi velocity of cold quasiparticles. 
However, for hot quasiparticles close to (ir,0), only the 
incoherent high energy feature around 200 meV is visible. 
The momentum dependence of this high energy feature, 
even though incoherent in its nature, is similar to the 
dispersion of a mean field SDW state: 



Et = \i^ + £k+Q ) ± \/ ( £k 2 £k+Q ) 2 + A§ DW , 

(40) 

where Ag DW = §A 2 , as can be obtained from the saddle 
point approximation of the Borel summed £ — > oo pertur- 
bation series (see appendix B). This provides an explicit 
demonstration that the high energy feature is indeed an 
incoherent precursor of an SDW state. The agreement 
between theory and experiment regarding the detailed 
momentum dependence of the high energy feature, is an 
important confirmation of the general concept of our ap- 
proach. 

While the overall position of the high energy feature 
(« 200 meV in the present case) depends on the value 
of t' , the general momentum dependence of these states 
remains robust against any reasonable variation of t' or 
the coupling constant g. We note that the experiment of 
Marshall et al. [[H] was performed in the strong pseudo- 
gap state. It is however natural to expect that the high 
energy feature remains unaffected by the opening of the 
low frequency leading edge gap; it will thus be the same 
in the weak and strong pseudogap state, and will be little 
affected by the superconducting transition. 

In ARPES experiments at half filling, it is found that 
the location of momentum states with half of the inten- 
sity of a completely occupied state, i.e., with — i, is 
nearly unchanged compared to the case at large doping. 
We show in Fig. |l5| our calculated results for the mo- 
mentum points with = i ; our results are quite similar 
for the physically relevant values 1 < £ < 4 and change 
only slightly for large correlation lengths. Note, that in 
the latter case nk varies only gradually. Even very close 
to (n, 0) it would be hard to determine experimentally 
whether nk=(7r,o) is larger or smaller than one half. Our 
results are therefore in agreement with the experimental 
situation; they demonstrate that there is a "memory" in 
the correlated system which, as far as the total charge of 
a given k-state is concerned, behaves quite similarly to 
the case without strong antiferromagnetic correlations. 

Finally, we address the question of why, for moder- 
ate values of the correlation length, we obtain such pro- 
nounced anomalies. In addition to £, the only length scale 
in the problem is the electronic length £o ss v/A 2w/g. 
It is natural to argue that once £ > £o some new behav- 
ior of the quasiparticles due to short range order might 
appear. Within standard weak coupling theories, £o is by 
construction a large quantity and the theory is trustwor- 
thy only for very large £. The summation of the entire 
perturbation series in our calculation however enables us 
to take account for the situation where £o can be of the 
order of a few lattice constants, i.e. for the intermediate 
coupling constant regime. That this qualitative argu- 
ment is also quantitatively correct, can be seen in Fig. 
where we show the £ dependence of the spectral density 
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at a hot spot for which e k = £ k +Q. For the above given 
set of parameters, £o ~ 2 and SDW precursors occur as 
soon as £ > £o- This is in striking agreement with, and 
provides a microscopic explanation for, the prediction by 
Barzykin and Pines that one finds £(T = T cr ) as 2 at the 
crossover temperature T cr , where the magnetic response 
changes character. 

We conclude that the quasiparticle excitations in the 
weak pseudogap regime are intermediate between a con- 
ventional system with a large Fermi surface and a spin 
density wave system with a small Fermi surface. The fact 
that both aspects are relevant explains the failure of any 
approach which concentrates on only one of these. 



VI. SPIN AND CHARGE VERTEX FUNCTIONS 

We turn now to the coupling of quasiparticles of the 
weak pseudogap state with the collective spin and charge 
degrees of freedom. This is of interest in its own right, 
and is of importance for an understanding of the charge 
and spin response functions discussed in II. The quan- 
tities which characterize the interaction of quasiparti- 
cles with the spin and charge degrees of freedom are 
the vertex functions r kk+q (w + iO + ,uj + v + i0 + ) and 
r k k+q (^ + i0 + , uj + v + i0 + ) . In order to have an idea of 
the behavior of these vertex functions we first consider 
their behavior analytically in the lowest nontrivial order 
of the perturbation series. Our subsequent numerical re- 
sults are obtained from the full solution of the problem. 
For the spin vertex, we find on using Eq. and Eq. || 
for S(q), that up to second order in g oc A: 



r k , k+q (w + * 0+ ' w + v + * 0+ ) = 1 



1 



3 uj — e k + iv/£ 

1 

w + v - e k+ Q + iv/£ 



(41) 



We are mostly interested in the vertex function for fre- 
quencies uj which correspond to the quasiparticle energies 
at the Fermi surface. For the case of an unchanged Fermi 
surface, the bare dispersion e k determines the quasipar- 
ticle energies at this Fermi surface. Once hole pockets 
are formed, these are given by the SDW energies E£ 
of Eq. Eo[ On evaluating Eq. [ll] at the SDW energies 
uj = S k of Eq. [lO] and for v = in the limit £ — > oo, 
we find, r s = |; the spin vertex is reduced. For the case 
of long range antiferromagnetic order, with only two spin 
degrees of freedom left, r s = vanishes, as was shown by 
Schrieffer Jij|. On the other hand, if one takes into ac- 
count that in the weak pseudogap regime the Fermi sur- 
face is basically unchanged and evaluates Eq. [ll] at small 
frequencies uj = v = for a hot spot with e k = £k+Q = 0, 
it follows that 



3v 2 



(42) 



i.e. the vertex is considerably enhanced. In the case of 
the charge vertex the prefactor | in Eq. [ll] has to be 
replaced by —1 and one finds r c = 4 if one considers 
£ — > oo and r c = 1 — ^-£ 2 in the case of an unchanged 
Fermi surface. These considerations demonstrate that 
which energies one considers and how the Fermi surface 
evolves is crucial for an understanding of the role of ver- 
tex corrections, i.e., enhancement vs. suppression. It 
also shows that only a careful and self consistent analy- 
sis can reveal in which way the renormalized charge and 
spin interactions vary. This we now do. 

In Fig. [l6]we show the spin vertex r k k+ q(^ + i0 + , uj + 
i0 + ) for hot and cold quasiparticles with momentum 
transfer Q and zero frequency transfer as function of 
energy and, for comparison, the corresponding spectral 
function, A k (w). As in the case for the spectral func- 
tion, the vertex function is strongly anisotropic; for cold 
quasiparticles vertex corrections are negligible, whereas 
the strong low frequency enhancement of T k k+ g (uj + 
i0 + ,uj + i0 + ) for hot quasiparticles demonstrates that 
despite their reduced low frequency spectral weight, hot 
quasiparticles interact strongly with the spin fluctua- 
tions. Thus, for physically reasonable ^-values, the low 
frequency vertex is not reduced but enhanced. This is im- 
portant and must be taken into account in constructing 
an effective theory for the low energy degrees of freedom 
of the strong pseudogap state. Even without detailed cal- 
culations, it is evident that in the spin fluctuation model 
the strong coupling nature of the low frequency degrees of 
freedom is crucial since it demonstrates that quasiparti- 
cles and spin fluctuations do not decouple. This is essen- 
tial to obtain a spin fluctuation induced superconducting 
state and suggests that new strong coupling phenomena 
are likely to occur once the temperature decreases and 
the system changes character due to the suppression of 
the quasi particle scattering rate. 

We further note that these results are frequency depen- 
dent. Thus when we consider the vertex at frequencies 
uj close to the high energy features, we find a moderate 
reduction, which has the same origin as the vanishing 
vertex of the long range ordered state p5| . 

The anisotropy of the spin vertex function can be seen 
in Fig. [l7|, where we plot the spin vertex r kk+ q(o; + 
i0 + , uj + i0 + ) for uj — along the Fermi surface as a func- 
tion of 4>f = axct&n(k y / k x ) (solid line). The rather sharp 
transition between hot and cold quasiparticle behavior is 
similar to that seen for the spectral density, in Fig. |lj. 

The corresponding behavior for the charge vertex func- 
tion is shown in Fig. [l8[ which shows that T k k+q (w + 
i0 + , uj + i0 + ) behaves in opposite fashion to the spin ver- 
tex. At low frequencies, the quasiparticles are almost 
completely decoupled from potential collective charge de- 
grees of freedom. This effect is strongest for hot quasipar- 
ticles, and occurs for momenta transfers q, around (0, 0) 
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as well as those close to Q. The formation of incoherent 
spin density wave precursors obviously leads to a decou- 
pling of the low energy quasiparticles from the charge 
degrees of freedom. In our theory, its origin is the domi- 
nant interaction of the quasiparticles with collective spin 
degrees of freedom. It is thus impossible that any kind of 
static or low frequency dynamical charge excitation can 
substantially affect the charge carrier dynamics of hot 
quasiparticle states in cuprates. Furthermore, this result 
explains another important puzzle of the cuprates: the 
irrelevance to transport phenomena in the normal state 
of the total electron phonon coupling, 

■< q Ph - * flfe" PMO) Il, k+q (0 + *0+, + i0+), (43) 

despite their pronounced ionic structure, which in fact 
suggests a strong bare interaction g£' q ph '(°) of charge 
carriers with the poorly screened lattice vibrations. 

These considerations show that an effective low energy 
theory of the weak pseudogap state must take the strong 
and anisotropic spin fermion interaction into account; it 
can safely neglect the coupling to phonons of hot quasi- 
particles as well as any other interaction they might have 
with charge excitations. This is an important new theo- 
retical constraint for the strong pseudogap state. 

VII. CONCLUSIONS 

We have used our solution of the quasistatic spin 
fermion model of antiferromagnetically correlated spin 
fluctuations to develop a new description of the interme- 
diate weak pseudogap state of underdoped cuprate su- 
perconductors. Based on the experimental observation 
that the characteristic energy scale of over-damped spin 
excitations is small compared to the temperature, once 
T lies between the two crossover scales T* and T cr , we 
conclude that the spin degrees of freedom behave quasi- 
statically i.e. the spin system is thermally excited and 
the scattering of quasiparticles with their own collective 
spin modes can be regarded as resulting from a static spin 
deformation potential, characterized only by the strength 
and spatial extent of these spin fluctuations. On neglect- 
ing the quantum dynamics of the spin modes, we were 
able to solve this spin fermion model by summing all di- 
agrams of the perturbation series for the single particle 
Green's function and the spin and charge vertex func- 
tions. This enabled us to directly investigate the spectral 
function, measured in angular resolved photo-emission 
experiments, the effective interactions of quasiparticles 
with spin as well as charge collective modes and the spin 
and charge response behavior. 

Our results demonstrate that for intermediate values 
of the antiferromagnetic correlation length and inter- 
mediate coupling constant g, one is dealing with the 
rich physics of a crossover between non-magnetic and 



spin density wave like behavior. While the Fermi sur- 
face of the quasiparticles remains unchanged their highly 
anisotropic effective interaction leads to two different 
classes of quasiparticles, hot and cold, with only the hot 
quasiparticles feeling the full strength of the antiferro- 
magnetic interaction. For the latter, a transfer of spectral 
weight to high energy features occurs. These high energy 
features are the incoherent precursors of a spin density 
wave state; their momentum dependence is in excellent 
agreement with corresponding high energy structures at 
around 200 meV seen in recent ARPES experiments. For 
low energies, hot quasiparticles have a reduced spectral 
weight, a weak pseudogap characteristic, and the coher- 
ent quasiparticle poles are completely over-damped due 
to the strong scattering rate. The drop in scattering rate 
found in the strong pseudogap state is not sufficient to 
make this quasiparticle pole visible. It is only below T c 
that the rate drops sufficiently that the pole becomes 
visible. This scenario explains the appearance of a sharp 
peak, invisible above T c , at low frequencies in the su- 
perconducting state of underdoped cuprates. The high 
energy features on the other hand are expected to be 
unchanged as temperature is lowered. 

Finally, we used our calculation of the irreducible ver- 
tex functions to investigate the effective interaction of the 
quasiparticles with spin and charge modes. We find that 
this effective interaction is likewise highly anisotropic; 
the low energy electron-spin fluctuation interaction is 
strongly enhanced whereas the coupling to charge de- 
grees of freedom is reduced. The enhancement of the spin 
vertex is essential for the development of a spin fluctua- 
tion induced superconducting state and is an indicator as 
well of anomalous behavior at lower temperatures. The 
reduction of the charge vertex causes a reduced electron- 
phonon coupling constant for hot quasiparticles as well as 
a decoupling of hot quasiparticles from potential charge 
collective modes. 

This scenario applies, as discussed, for temperatures 
between and T cr . Below X 1 *, the characteristic fre- 
quency of the spin system increases for decreasing tem- 
perature, making it impossible to consider the system 
as exhibiting quasistatic behavior. Thus, here we ex- 
pect that the quantum nature of the spin degrees of free- 
dom becomes essential; it brings about a strongly re- 
duced phase space for the inelastic scattering of quasi- 
particles and spin fluctuations, and causes a sudden drop 
in the corresponding scattering rate of the system. Nev- 
ertheless, since we expect the high energy features to re- 
main unchanged, the spectral weight of hot quasiparticles 
is strongly reduced while their interaction is enhanced. 
This strong coupling behavior can bring about precursor 
to superconducting state behavior, as recently discussed 
by Chubukov |46[ , who found that at zero temperature, 
the quantum behavior of the spin modes causes a lead- 
ing edge gap of the spectral function within a magnetic 
scenario for cuprates. It has recently been shown that 
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this leading edge gap will be filled with states once the 
temperature is raised, and that it basically disappears 
for T w 0.6w sf £ 1/l/ , with i/«2§. Above this tempera- 
ture, the system starts to behave quasistatically and the 
superconducting precursors are suppressed. Within the 
same scenario, one thus finds that the classical behav- 
ior of spin fluctuations favors magnetic precursors and 
causes an increase of the low energy interaction, which 
leads for lower temperatures to superconducting precur- 
sors due to the quantum behavior of the spin fluctua- 
tions. This provides perhaps the first consistent picture 
of the microscopic origin of the two pseudogap regimes 
and their crossover temperatures: The strong pseudogap 
state and the leading edge gap are the quantum manifes- 
tation of strong antiferromagnetic correlations whereas 
the spin density wave precursors are its classical mani- 
festation. 

What are the consequences of this new microscopic sce- 
nario for the crossover behavior of underdoped cuprates? 
Our theory predicts that upon measuring the uniform 
susceptibility Xo(T) and the spectral function for an 
slightly underdoped or optimally doped system, (for 
which T cr is not too high), it will turn out that above T cr , 
where Xo(T) is maximal, the high energy features will dis- 
appear. We also expect important insights into the role 
of impurities and high magnetic fields. Non-magnetic 
impurities should affect the weak pseudogap state only 
slightly. Their strongest effect will occur slightly below 
T cr where impurities cause £ to decrease to £i mp . Pro- 
vided £i mp < £ » 2 < £, an impurity driven transi- 
tion out of the weak pseudogap state occurs. T* will 
be considerably more sensitive to impurities because for 
particle-particle excitations with a tendency to d-wave 
pairing, non-magnetic impurities act destructively and 
the strong pseudogap state may even completely disap- 
pear. 

High magnetic fields provide another indicator of the 
differences between strong and weak pseudogap behav- 
ior. In the weak pseudogap regime no sensible effects 
for all achievable field strengths will occur, because the 
relevant coupling of the magnetic field is the Zeeman in- 
teraction with the spins, oc H • s, which has to compete 
with the much stronger short range correlations of these 
spins. In the strong pseudogap state, we expect that the 
dephasing due to the minimal coupling p — > p — |A will 
strongly affect the behavior in the pairing channel, caus- 
ing a suppression of strong pseudogap behavior. From 
this perspective, it immediately follows that the trans- 
port experiments by Ando et al. |49[ ] in a pulsed high 
magnetic field demonstrate that weak pseudogap behav- 
ior is prolonged to lower temperatures. It is tempting to 
speculate that for very strong magnetic fields the weak 
pseudogap regime crosses directly over into an insulating 
one. 

In general , for decreasing temperature the weak pseu- 
dogap regime can therefore cross over to a strong pseu- 



dogap state or to an insulator and, under certain cir- 
cumstances, directly into the superconducting state. The 
latter may occur for systems with very large incommen- 
suration, which tend to reduce the scattering rate of a 
quasistatic spin system. [ p0[ ; it is likely of relevance for 
La 2 _ a; Sr a ;Cu04. 

Our results for the low frequency spin dynamics and 
optical response as well as the Raman intensity in the 
Bi 5 channel will be compared with experiment in a sub- 
sequent paper. There we also demonstrate that we can 
explain the generic changes of the low frequency magnetic 
response at the upper crossover temperature T cr pl| . 
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APPENDIX A: SINGLE PARTICLE GREEN'S 
FUNCTION AND DIAGRAMMATIC RULES 

In this appendix we derive the diagrammatic rules of 
the spin fermion model under circumstances that the 
interaction between collective spin modes is neglected. 
Note, these diagrammatic rules do not rely on the quasi- 
static approximation, but are completely general for the 
spin fermion model. 

Using the generating functional 

W[r], rf\ = \ J VdVc exp {-{S + c+tj + 77+c)} (Al) 

with partition function Z = j'Dc^'Dce~ , effective 
action of Eq. fj] and short hand notation c*rj = 



Jo c ker r 7k,<T! the single particle Green's function 



can be obtained via a functional derivative: 



HL( t )<W(t') 



(A2) 



77=77^=0 



As usual c, etc. are Grassman variables. In order 
to perform this functional derivative it is convenient to 
introduce a collective bosonic spin 1 field S q by adding 
an irrelevant Gaussian term to the action S — > S + S , 
where 
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S (S) = \[dr j* dr' X^(r - r') S q (r) • S_ q (r') , 

(A3) 

to integrate with respect to this spin field and to divide 
by the corresponding partition function Zb of this ideal 
Bose gas. Finally, we shift (after Fourier transformation 
from time to frequency) the variable of integration as: 



Sq(«W„) -> Sq(iu}„ 



2g_ 

V3 



X q («w n )s q (iw„) , (A4) 



leading to the effective action of the spin fermion model: 



. c L G ok c k<' 



L dr ^ 

J ° k , <j 

^|^rfr^s q (r).S_ q (r) 



(A5) 



After this Hubbard Stratonovich transformation, we can 
integrate out the fermions, yielding 



W[n,r, ] ] 



1 



zz B 



PSexp 



(A6) 



with the action of the collective spin degrees of freedom 
S C (S) = -trln (-G(S)- 1 ) + S (S) . (A7) 
Here, we have introduced the matrix Green's function 

--S_S k _ k /(T)*(r-T').<7, (A8) 

which describes the propagation of an electron for a given 
configuration S of the spin field. Performing the above 
functional derivative with respect to r\ and rf gives finally 



T ) = — \ GkMaa(T,T'\S) 



Gk.(r 

exp{trln(-G(S)" 1 )}J) , 

where the average 

(•••>»= ^ / I?S-..exp{-5 } 



(A9) 



(A10) 



is performed with respect to the free collective action of 
Eq. A3. This is a standard exact reformulation of Eq. |?] 
of collective spin fields which has the appeal that explicit 
fermion degrees no longer occur and that Wick's theo- 
rem of the Bose field S can be used to evaluate the single 
particle Green's function. 

Since we do not expect the interaction of the spin 
modes to be relevant, we neglect nonlinear (higher order 



in S than quadratic) terms of the spin field, assuming 
that no modifications due to spin fluctuation-spin fluc- 
tuation interactions occur beyond those already included 
in Xq(u)). Using renormalization group arguments, it has 
been shown by Millis |5l| that indeed the system is char- 
acterized by a Gaussian fixed point as far as the low fre- 
quency spin dynamics is concerned. The mathematical 
consequence of these two assumptions is that we can use 
the approximation 



exp{trln(-G(S)- 1 )} a Z. 



(All) 



Contributions of second order in S can be ignored, since 
they would renormalize Xq(^), which is assumed to be 
the experimentally determined, i.e. fully renormalized 
susceptibility. Eq. All leads to a considerably simpli- 



fied expression for the Greens function: 



G kCT (r-r') = (G k , kffff (r,r'|S) 



(A12) 



Consequently, the diagrammatic series for the deter- 
mination of the single particle Green's function reduces 
to that of a single particle problem with time dependent 
spin " impurities" . Inversion of the Greens function ma- 
trix of Eq. A8 in spin space yields for the cr-spin matrix- 
element: 



G<j(S) — ^G Q 1 — ~J=S Z 



V3 



(A13) 



with S a = S if er = ±1. Here, we still have to take the 
matrix nature of G D and S (in momentum and frequency 
space) into account since they are diagonal in different 
representations (G in momentum and frequency space, 
S in coordinate and time space). Considering the limit- 
ing case of only one longitudinal spin mode generated by 
S z , gives after averaging: 



(G ; (S)) -G ^ 



N=0 



2\ N 



(S GqS Go) 



N 



(A14) 



where we used the fact that odd orders in S z vanish 
without global symmetry breaking. Using this represen- 
tation of the Greens function, the diagrammatic rules 
which correspond to the averaging with respect to S z 
follow straightforwardly. The averages can be evalu- 
ated via contractions based on Wick's theorem using 
{S^{t)S z _^{t')) = x q (r - t'). This leads to the dia- 
grams of a theory of fermions interacting with a scalar 
time dependent field. Here, the topology of each diagram 
is identical to the topology of the contraction symbols 
which occur by applying Wick's theorem. Alternatively, 
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one can also consider the case of two transverse modes 
leading to 



(G t (S)) = G J2 



N=0 



'S-G S+G ) N ) (A15) 



which is identical to a theory of fermions interacting with 
a "charged" time dependent field. Again Wick's decom- 
positions using (Sq (t)SZ ci (t'))o — 2x q (T — t') can be 
performed. 

The situation becomes more complicated if one con- 
siders simultaneously lon gitud inal and transverse modes. 
Here, it follows from Eq. A13: 



G & (S) = GoAS') £ ( 

N=0 V 6 7 



N 



x (S- G ,- a (S z )S+ GoAS z )) 



N 



where 



G ,a{S z ) = G £ ( a ) N ( -^S Z G 



/V 



(A16) 



(A17) 



Eq. |A16| and |A17| mean that for any transverse (spin flip) 
scattering event all possible longitudinal (spin conserv- 
ing) processes occur. The complication of this inserted 
partial summation is the occurrence of the sign factor 
(cr) of down spins. The (— 1)^ factor occurs, if there 
are contractions out of longitudinal processes for an odd 
number of S z fields of a a =1= —1 Green's function 
G 0t i(S z ). If a S z is paired with another S z of the same 
Go,i(S z ), they occur in an even number without modi- 
fying the sign. The (— 1) N enters only if a longitudinal 
field S z is contracted with one which refers to a Greens 
function G ^(S Z ). In order to reach G 0} ^(S Z ), the corre- 
sponding line of the S z contraction has to cross an odd 
number of transverse contractions. From these consider- 
ations, the following diagrammatic rules for an arbitrary 
diagram of order 27V result: 

• Draw 2 TV + 1 solid lines with 27V vertices which can 
be 2L spin conserving, TV — L spin lowering and 
TV — L spin raising vertices , referring to vertices of 
longitudinal processes (S z ), leaving transverse pro- 
cesses (S~) and entering transverse processes (S + ), 
respectively. For the transverse vertices one has to 
ensure that two subsequent spin raising (lowering) 
vertices are separated by one spin lowering (rais- 
ing) and an arbitrary number of spin conserving 
vertices. 

• Connect the spin lowering vertices pairwise with 
spin raising ones by a wiggly line 

• Connect the spin conserving vertices pairwise with 
dashed lines 



• Insert for any solid line a Green's function 
G k(icu n ), for any vertex, g/y/3, for any wiggly or 
dashed line x q (iw„) and take momentum and en- 
ergy conservation at the vertices into account. 

• Multiply with a prefactor 2 N ~ L which accounts for 
the two transverse spin degrees of freedom (TV — L 
is the number of wiggly lines). 

• Multiply with a prefactor (— l) c in front of the dia- 
gram, where C is the number of crossings of dashed 
and wiggly lines. 

Finally one has to sum over all possible diagrams gener- 
ated by this procedure. In Fig. |l^ we plot all diagrams for 
the self energy up to order g 4 , including signs and multi- 
plicities. The occurrence of the additional crossing sign 
which results from the interference of longitudinal and 
transverse modes certainly complicates the situation. If 
there were only one or two components of the SU(2) spin 
vector, the diagrammatic rules reduce to the special case 
of only longitudinal or transverse interaction lines. Here, 
the problem is identical to that of uncharged or charged 
bosons, respectively. Only the simultaneous considera- 
tion of both phenomena, which is necessary to preserve 
spin rotation invariance, leads to the prefactor (—1) and 
reflects the fact that the collective mode of the system is 
a spin fluctuation. 

Identical diagrammatic rules can be derived for the 
spin fluctuation vertex function. In Fig. ^ we show, as 
an example, the spin flip vertex T t as well as the spin 
conserving vertex T\. Even though our approach is not 
constructed to be manifestly spin rotation invariant, this 
symmetry must of course be fulfilled once the transverse 
and longitudinal spin susceptibilities, represented by the 
wiggly and dashed lines of the above diagrammatic rules, 
are the same. As can be seen from the lowest order vertex 
corrections, the above derived diagrammatic rules guar- 
antee indeed that T\ = T t = r s , as expected. It can 
be shown order by order in the coupling constant g that 
our procedure guarantees spin rotation invariance, as is 
essential for a system with isotropic spin fluctuations. 



APPENDIX B: DETERMINATION OF THE 
MULTIPLICITY OF DIAGRAMS 

In this appendix we calculate the multiplicity of iden- 
tical diagrams for a given order of the perturbation se- 
ries. This will be done in two steps. First, we solve the 
problem in the limit £ — * oo; second, we use the general 
expression of Eq. valid for an arbitrary diagram and 
finite £ and determine the missing ^-independent multi- 
plicity factors from the £ — > oo solution. For the special 
cases of only longitudinal or transverse modes, our solu- 
tion is the same as Sadovskii's |29). It is important to 
notice, that without the result of Eq. it would not be 
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possible to determine uniquely the diagram multiplicity 
from the infinite £ limit. 



and 2, the analytical inversion of the Dyson equation is 
evident, for Af = 3 it follows: 



a. Solution for £ — > oo 

The limit £ — * oo is not free of complications: First, 
we expect that in this limit the longitudinal and trans- 
verse spin degrees behave differently. We can ignore this 
problem here because we are only interested in the spin 
rotation invariant situation for finite £ and use the limit 
only for the mathematical purpose of determining dia- 
gram multiplicities. Thus we assume spin rotation in- 
variance also for £ — > oo. Second, the local moment of 
the susceptibility in Eq. diverges in the static limit loga- 
rithmically for £ — ► oo. This problem can also be avoided, 
because the use of Eq. |l6| avoids this divergence, but does 
not change the multiplicity of the diagrams. Third, the 
rather straightforward result that for a given order in g 
each diagram for £ — > oo is besides sign and multiplicity 
identical, leads to the following perturbation expansion 
of the Green's function: 



oo 



GkM = G ok > j (2n 



n=0 



'A 2 

M"[- 



k > 



(Bl) 



where = G kG k+Q, which is in fact a divergent se- 
ries. One can then obtain a convergent result using Borel 
summation of this series. This is however not the most 
transparent way to solve this problem; we choose an alter- 
native approach, using the path integral representation 
of the Green's function derived in appendix A, which of 
course gives the same result. 

In the limit of infinite antiferromagnetic correlation 
length, the only relevant spin configuration is: 



k-k' 



S5 



k,k'+Q 



(B2) 



and the path integral of Eq. AlC simplifies considerably 
to 



- u - 

2A 2 



(B3) 



Here, the spin rotation invariant case of present physical 
interest is Af = 3. If Af = 1 or 2 the system consists only 
of longitudinal or transverse modes, respectively. In the 
following we solve the problem for all three situations. 
In doing so, we have to replace in all expressions g 2 /3 by 
g 2 /Af. For example one has to generalize the expression 



* 2 = ^ 2 



(B4) 



for the gap energy etc. Here, the partition sum of the 
bosons is given by Z B = ±r Ar r(Af/2)(^-) Ar / 2 with 
n = 2, r 2 = 2n and r 3 = An, respectively. For Af = 1 



g ct (s; 



k,k' 



1 



G k 



-<5k,k' + er- 



1 



\/3 



?k,k'+Q ■ 



3 "k 3 

The second, non-spin rotation invariant, term vanishes 



after averaging, Eq. (B3), and we can finally write: 



GkM 



,W-£k" 



(gg) 2 /3 



(B5) 



It follows that the full Green's function, obtained in the 
limit £ — > oo, is an averaged second order Green's func- 
tion with fluctuating SDW gap A = y/ (gS) 2 /3. For ar- 
bitrary Af, the same result occurs if one replaces the gap 
by A = J (gS) 2 /Af. Performing the angular integration 
of the vector S, the integral of Eq. B3 can be written as: 



GkM = 




(B6) 



The remaining one dimensional integral demonstrates 
the different behavior for different Af. For Af = 1, 
the distribution function of the SDW-gap is centered 
around zero whereas for Af = 2 and 3 it has a maxi- 
mum for finite A. Performing the saddle point approx- 
imation for Af = 2 or 3 yields SDW-like solutions with 

reduced gap A c = yJ^A w 0.7071A for Af = 2 and 

A = y/|A « 0.8165A for Af = 3. However, even for 
Af = 1, the contribution of the tails of the distribution 
function changes the behavior qualitatively compared to 
the saddle point approximation and solutions similar to 
an SDW state occur. It is interesting to note Eq. B6 is 



the Borel integral representation of the formally diver- 
gent perturbation series in the limit £ — > oo. With the 
path integral approach, we did not encounter this diver- 
gence, thus demonstrating that it was, in fact, a spurious 
one. 

The integral with respect to the fluctuating SDW gap 
can be evaluated using the integral representation of the 
incomplete Gamma function 



r(i-v)y z + t 



■dt. 



(B7) 



with t 



= AAA 2 /(2A 2 ) and ip = 1 - Af/2 and z = 
£k)(uj — £k+q)/A 2 . Using the continued fraction 
representation of T(ip, z) [ p2[ , we obtain the following re- 
sult for the single particle Greens function: 
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W — £k+Q 



« 2 A 2 



- £ k 



was approximated by Eq. 21 The decoupling of the mo- 
mentum integrals in the proper time representation and 
the diagram multiplicities, i.e. the Kj, do not depend 
on the actual choice of ■f/'k+Q^): and the only difference 



w — £k+Q - • • • 
(B8) 



compared to Eq. B1C is the function g^\u>), which in 
terms of -0k+Q(£) can be expressed as 



poo 

with k 3 = j/N if j even and Kj = (j +J\f- if j odd. g U) ^ = -i dt e ^~^+^) </, k+jQ (t) 3 . (B12) 

For the special cases of M = 1 and 2, this is identical to ^ Jo 



For the special cases of N = 1 and 2, this is identical to 
Sadovskii's result |29| ]. Furthermore, we can obtain an 
analytical expression for the single particle Greens func- 
tion for the case N = 3 which corresponds to the spin 
fermion model. 



On using the result of Eq. [2C] for -0k+Q(i), together with 
Eq. BIO, the solution of the Green's function of the spin 



fermion model is that given in Eq. 3C 



b. Generalization to the case of finite £ 



Using, for the moment, the approximation of Eq. |T^, 
the solution for finite antifcrromagnetic correlation 
length can be inferred from Eq. B8 and Eq. valid 
for an arbitrary diagram and finite £. From Eq. ^9] we 
know that, compared to the limit £ — > oo, the only way 
the correlation length enters the problem is via: 



e k+jQ 



(B9) 



where j refers to the order of the continued fraction with 
nominator KjA 2 . On the other hand, in Eq. [2^, the in- 
teger number rij has a specific diagrammatic meaning. 
Since we can generate an arbitrary diagram by expand- 
ing the continued fraction of Eq. B£ with respect to A 2 , 



we can perform the replacement of Eq. BE within the 
continued fraction representation Eq. 



In order to fix 

the not yet determined integer number rij of the j-th in- 
sertion of the continued fraction, which is independent of 
the correlation length, we use the fact that for a given 
order A 2N of the perturbation theory, only one sequence 
{rij} occurs with rij = N . This refers to diagrams which 
are identical to the rainbow diagram of order 2N . The 
only way that this can be fulfilled is rij = j. Thus we 
obtain the solution of the spin fermion model for finite 
correlation length based on the approximation of Eq. [If]: 



iA 2 G^ +1 \uj). (BIO) 



with 



<?k j) M- 



W - Ek+jQ + i^f^ 



(Bll) 



which generates the continued fraction representation of 
the single particle Green's function, similar to the one 
dimensional case. 

The general solution, independent of Eq. |l^, follows 
from the fact that the only difference from the strict two 
dimensional case is the function ^k+Q (t) of Eq. n9l which 
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FIG. 1. Schematic phase diagram of underdoped cuprates. 
Note the presence of two different crossover temperatures: 
T CI , which characterizes the onset of sizable antiferromag- 
netic correlations; and T* , which signals the onset of a con- 
siderable loss of low energy spectral weight in the quasipar- 
ticle spectrum leading to a minimum of the characteristic 
spin-fluctuation energy u> a f ■ The region between T cr and T* 
is the weak pseudogap regime discussed in this paper. 
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FIG. 2. ARPES spectra from Ref.[30] for momenta k at 
and close to (tt, 0), for two different doping concentrations. The 
T c = 78 K sample is slightly overdoped whereas the T c = 88 K 
sample is underdoped. 
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FIG. 3. The position of local maxima of the spectral func- 
tion along the high symmetry lines of the Brillouin zone is 
shown for an overdoped and underdoped system (data from 
Ref. [31]). 
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FIG. 4. Fermi surface for an overdoped and underdoped 
system, obtained from k-points where local maxima of the 
spectral function cross the Fermi energy (data from Ref. [31]). 
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FIG. 5. Spectral function of a hot quasiparticle above and 
below the superconducting transition temperature (data from 
Ref. [22]). 
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FIG. 6. A typical bare Fermi surface in the first quarter of 
the BZ, closed around the momentum point (tt, 7r), for differ- 
ent AF correlation lengths. The thick sections characterize 
the hot parts of the Fermi surface. 




FIG. 7. Illustration of the sequence {rij} for a self energy 
diagram of order g 8 . rij is the number of spin fluctuation lines 
above the j-th fermionic Green's function. 
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FIG. 8. Self energy diagrams of order g 4 and g 8 , which are, 
within the quasistatic approximation, identical apart from 
multiplicity and sign, because the number of spin fluctuation 
lines (shown by wiggly lines for transverse and dashed lines 
for longitudinal spin excitations) on top of a given electron 
propagator (solid line) are the same. 




FIG. 9. Diagrammatic illustration of the recursion relation 
of the vertex function, which is similar to the ladder approx- 
imation of the irreducible vertex. All non-ladder diagrams 
are taken into account by the corresponding weighting fac- 
tors ri+i. 
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FIG. 10. The spectral density multiplied with Fermi func- 
tion on the Fermi surface for £ = 3. The distinct behavior of 
hot and cold quasiparticles is visible. 



(0,0) (O,,!) 

FIG. 11. Fermi surface (solid line), in comparison with 
the bare Fermi surface (dotted line) and the visible part of the 
Fermi surface (thick solid line) , i.e only for momenta where a 
maximum of the spectral density crosses the Fermi energy if 
k crosses the Fermi surface. 
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FIG. 12. The spectral density ^4k(^ = 0) along the Fermi 
surface as function of 0f = a,rctan(k y / k x ) . 




FIG. 13. the momentum dependence of local maxima of 
the spectral density as function of £ and hole doping con- 
centration nh is compared with the experiments of Ref. [25] 
for Bi2Sr2Cai_ I Dy I Cu208+a with x — 0.1 (triangles) and 
x — 0.175 (diamonds). Only maxima with relative spectral 
weight > 10% are shown. 
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FIG. 14. £ dependence of the spectral density at the hot 
spot. The appearance of SDW precursors for £ > 2 can be 
seen. For smaller values of £ the system behaves convention- 
ally. 




FIG. 15. k -points with rik = \ for different correlation 
lengths in comparison with the results for an uncorrelated 
Fermi system. 
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FIG. 16. The spectral density A^(ui) and spin vertex func- 
tion r^k+Q^ + i0 + , lo + i0 + ) for a hot quasiparticle with mo- 
mentum transfer Q and zero frequency transfer as function of 
energy are shown for two correlation lengths. 
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FIG. 18. Charge vertex T^ k+q (u + 10 + , lo + i0 + ) for lo = 
as a function of <j)F = axcta,n(k y /k x ). 




FIG. 19. Self energy diagrams of the spin fermion model 
up to order g 4 , including signs and multiplicities. The wig- 
gly (dashed) lines correspond to spin flip (spin conserving) 
processes. Each crossing of a wiggly and dashed lines causes 
a prefactor —1 and each wiggly line an additional factor of 
2. The solid line corresponds to the bare fermion propagator 
and each vertex to the coupling constant g/3 1 ^ 2 . 
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FIG. 20. Longitudinal (spin conserving) and transverse 
(spin flip) vertices up to g 2 . Note, that the derived diagram- 
matic rules guarantee, as expected, spin rotation invariance: 

r, = iv 
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